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Vortex Motion in certain Triangles. 

By A. B. H. Love, B. A., Fellow of St. JohrJs College, Cambridge. 



1. In the Proceedings of the London Mathematical Society, Vol. XII, Routh 
has explained a method whereby the motion of a fluid in two dimensions due to 
sources or vortices situated within a given region can be inferred from the motion 
due to sources or vortices of equal strength placed at the corresponding points of 
another region, whenever, by means of a transformation by conjugate functions, 
we can obtain a correlation* ("conforme Abbildung") of the first region with the 
second in such a way that elementary portions of the two regions are similar. 
The conditions for such correlation may be expressed as follows : Let Z be a 
function of a complex variable z, and suppose Z and z represented in the usual 
manner by points on two planes, then the boundary of any region in the z plane 
will be transformed into a certain curve in the Z plane. The function Z is to be 
chosen to be finite, continuous, and one-valued within the region in the z plane, 
and points on the boundary in the z plane are to correspond to points on the 
boundary in the Z plane in such a way that one point on either boundary corres- 
ponds to one point on the other. 

If <3> and I 1 are the velocity- and stream-functions of any fluid motion in the 
Z plane,- then <E> + <,*P is a function of Z, and Routh shows that the velocity- and 
stream-functions in the corresponding motion within the correlated region of the 
Z plane can be inferred by simply substituting for Z in <E> + 0$ its value expressed 
as a function of z. 

In Art. 21 of Riemann's "Inaugural Dissertation" (Ges. Werke, p. 40) it is 
stated that the problem to find a correlation of two simply-connected plane regions 



♦For the theory of the correlation of two plane regions the reader may consult Darboux, " Theorie 
Generate des Surfaces," ch. IV, or Holzmiiller, "Einfuhrung in die Theorie der isogonalen Verwandt- 
schaften und der conformen Abbildungen. " 
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has always one and only one solution,* and thus in particular it is always 
theoretically possible to correlate any given simply-connected plane region with 
that part of a plane which lies on one side of the axis of real quantities. In a 
memoir in Borchardt's Journal (Bd. 70), Schwarz has discussed the solution of the 
problem when the region to be correlated with the half plane is bounded by a 
polygon, and he has shown that when the polygon is a triangle in the z plane, the 
function Z is in general a transcendent having an infinite number of values at 
every point of the z plane, of which the value along any given branch of the 
function is, however, singly determinate (eindeutig) within the triangle. But there 
are certain exceptional cases in which Z is an algebraic function of an elliptic 
function ^ of z, pz being one-valued at every point of the z plane. These excep- 
tional cases are those of a triangle whose angles are 

fit 71 7t\ (it 71 7t\ /27t Tt 71 \ / K 7t 7t\ 

\T ' T ' YJ ' VT ' 3" ' V/ ' Vs ' T ' TV ' or VT ' T ' Tv 1 ' 

The integration of the differential equation on which the correlation depends 
had been considered by Briot and Bouquet in their "Theorie des Fonctions 
doublement periodiques."f 

In what follows I propose to work out in detail these cases, and to apply the 
solution to the problem of the motion of a single vortex within the triangle. 
The solution for any number of vortices can be found by summation, while that 
for sources is obtained by interchanging <J> and I s . 

2. Schwarz shows that if a, b, c are points on the real axis in the Z plane, 
then a triangle whose angles are arc, fire, yit in the z plane can be correlated 
with the half plane above the real axis in the Z plane by means of the equation 



d ( 
dZ\ 



l ° g dz) Z—a + Z—b + Z—c 

and the corners of the triangle will correspond to the points Z= a, b, c. 

This appears because (1) the function z of Z is plainly in any correlation of 
the form C x z-\- 2 , where G lt G 3 are arbitrary constants, and z is one determina- 
tion of the function ; ( 2) the transformation 

d /, _dz\ a — 1 



log 



dZ\ & dZJ Z—a 



*A11 solutions are counted as one solution in which finite, triangles remain similar. 
t See the same author's "Theorie des Fonctions Elliptiques," pp. 388 seq. 
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effects a correlation of the half plane Z (imaginary part positive) with the space 
bounded by two lines in the z plane inclined to each other at angle an, and 
meeting in the point corresponding to Z=- a. 

It follows that the desired correlation is effected by means of the equation 

c lS + <7 2 =J % \z— a y-\z— hf-\z— c y- i dz, (A) 

in which the points a, b, c are taken to lie in a definite order on the real axis, 
viz. it A, B, G be the order of the corners of the triangle in the z plane at 
which the angles are an, fin, yn when its contour is described in the positive 
sense, i. e. with its area on the left, then a, b, c is to be the order of the points 
in the Z plane which correspond to the corners, the real axis being described in 
the positive sense. 

Case I. — The equilateral triangle. 
3. In this case the equation (A) becomes 

a% + G % = r n — d ^~~* a • C 1 ) 

1 T 2 J {Z—a)\Z—b)\Z—6f W 

We may take a = , 6=1, c = — 1 ; then, supposing Z= — ■=■ , we have 



/z' 



dZ' 



Wrlte l-Z' 2 = ^ = l--i-, (2) 

thus 



(Z' a — l)t 
1 



C x z+C % --—j V1 l a;8 , 
Hence, adjusting constants, we may write 
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(3) 



(4) 
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The correlation is effected by means of the equation 

1 3 

giving z% _ 1 — 1 



\—p\ (pz — e0( pz — %)(pz — <%) ' 
thus iZ _ 1 , ■> 

4. The elliptic functions given by (4) have been discussed by G-reenhill.* 
We have to notice that the discriminant is negative and equal to — 432. The 
real half-period o 2 is given by the equation 

/ *» ax _ 1 \67 \~2~) 
ViV4(^-l)~ 6 r /_2\ ' 



dx 
G) 2 



the corresponding imaginary half-period a 2 = ia 2 \/3. The periods a lt a 3 corres- 
ponding to e 1 , e 3 are % = — (o 2 — Og') , <-> 3 = — (o 2 + o0. 

We can also express in a simple form all the roots of the function pz. We 
know that for pure imaginary values of z, pz passes through all real values 
between e 2 (= 1) and — oo, and takes the same value at points whose z differ by 
multiples of 2o) 2 . Hence there is a root for some value of z a pure imaginary 
between and 2o 2 . Now we may show that if pz = 0, then p(Sz) = co , this 
follows at once from the formula 

p n z{2ph — 40p»z — 16) 



pBz — pz — 



3(£>*2 — 4f*) 



2 
by putting pz = . Thus it appears that -r-Wj is a root of pz which is pure 

o 

4 
imaginary and lies between and 2o 2 ', the other root -^-o 2 may be regarded as 

2 
— o 2 -f- 26) 2 , so that it differs from the former by a period and a change of 

o 

sign. 

5. Returning to equation (5), we see that p s z = gives the values of z which 
correspond to Z = ± 1 , and pz = oo corresponds to Z = . We may choose 

*"The Trajectory for the Cubic Law of Resistance," Proc. Royal Artillery Institution, 1886. 
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the pole z for Z=0 to be z = 0, and the root z for Z = — 1 to be z = — «;, 



2 
~3~ 



then since — o^ = (c^ + — af) — 2o 3 , we have p fa 2 + -=- of) = 0, and 

s = &) 3 + -5- «2 may be taken to correspond to Z = 1 . 

The function if= 2t/p'a is finite, continuous, and one-valued within the 
equilateral triangle in the z plane whose corners have coordinates 

( ' )'( '^°»)' (" 2 '^3" 2 )' 

and the points on the real axis in the Z plane correspond to points on the sides 
of the triangle, the points Z— 0, 1, — 1 being those which correspond to the 
corners. The height of the triangle is the real half-period a> 2 . 




-X 



-1 



Z 



The figure represents the boundaries of the two correlated regions, the points 
0, 1,-1 correspond to the angular points of the' triangle and these are marked 
with the same numbers in both figures. 

6. Now let z be the point at which is a vortex of strength m, and let Z be 
the corresponding point of the Z plane, the imaginary part of Z being supposed 
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positive. Then if Z' a , z' be the conjugate imaginaries of Z a , z , the fluid motion 
is given by 

T 2«- *Z—Zl 2n & p> z — p% K ' 

The position of z„ is given by p\ = 1 — -^-„, and if z be one of the roots of 

this equation, there will be other roots which do not differ from z by multiples 
of the periods. There are in fact two such roots %, z % for which pz 1 = e 1 pz 0) 
pz % = e 3 pz , and then z 1 =e 1 z and 2% = e 3 z . 
This appears from the equation 



OT^-d^kw-u 



e a being a cube root of unity, from which we deduce 

petf — eipz) 
so j3e s 2 = e 8 j3aj 



(7)* 



The roots z , %, z 2 of the equation p 3 z = 1 — -=5. satisfy the relation z + s x 

+ 2a = 0. If gj be the conjugate imaginary to %, then the conjugate imaginaries 
z{, zl to %, z % are e 3 2o, e^. 

Now consider the function 



ft (z)= 



1 

1 ^ p'^ =(jWj— ^Xp'z — p'sso), (8) 

1 



pz 


p's 


JWl 


f'si 


F2» 


§»'«» 



where p'^ = p% = pV 
We know that 



/ \ „ S fa — z) (5 (z 2 — z) (3 (z + z,+ Zo) <3 (z 2 — Zt)+ . . 

♦» <*> = 2 6wU 8 T • ( 9 ) 



Hence, remembering that 63 is an odd function, and that 

_ _ (3 (z 2 -j- z x ) (5 (z 2 — zj) 



ft=- WT ^^ . do) 



♦Greenhill, loc. cit. Art. 6. 

t Halphen, "Traite des Fonctions Elliptiques," p. S19. 

22 
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we obtain 

9 t z _ p , - 2 <S(z—z )6(z-z 1 )6(z—z 2 ) ^ 

&Z(3Z^3Z 1 <3Z % ' ^ ' 

so that the fluid motion is given by 

*^ tV 2n i0 ^ <S(z-z' )6(z-z[)<o{z—z',)' ^ 12) 

7. Greenhill* has considered this case of fluid motion proceeding by the 
method of images. He observes that the images are 

and the similar sets grouped round centres of hexagons whose coordinates are 
2m.au, 2»?(v/3; (2m + l)o 2 , (2m + 1)g) 2 V3. As 'the vectors of all these are 
included in the formula 2mo> 1 + 2m / co a , the solution found by summation will be 
that in (12), since 

<oz = znl(l -~)e» +ir ™ 2 , ' 

where w has all the values taken by the quantity 2?na 1 + 2m'a s , when m, m' are 
integers varying from to ± oo , except the value 0. In Greenhill's solution, 
which differs from (12) in form, the 6 functions are formed with a 2 and ol as 
fundamental half-periods, and, since u z + u' % is a period, the S 2 so formed is 
included in the 6 formed with Wj , w 3 as fundamental half-periods. 

Case II. — The right-angled triangle containing an angle of 60°. 

8. In this case equation (A) becomes 

z 

{Z—a)\Z—b)\Z—cf 
~Z 



OlZ+Gs—f— -y= TTT— -,. (13) 

J (Z — a)\Z — by{Z — c)<? 
Take a = , b = 1 , c = — 1, and suppose z = -=- , then 



Ot+^-fxz* 



dZ> 



(l— zy(ij-z')V 

Write ^Zl.^a-iz^ , , 

*" Applications of Weierstrass's Elliptic Functions," Proc. Lond. Math. Soc, Vol. XVIII, p. 377. 
t Schwarz, " Formeln und Lehrsatze zurn Gebrauche der elliptischen Functionen," p. 5. 
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then 



„ - p* dx 



dx 



1) 



(15) 



Hence, adjusting constants, 

x = pz, (16) 

where the invariants and periods of the elliptic functions are the same as in 
Case I. 

The correlation is effected by means of the equation 

1 — Z _ 3 

i + z~ pz ' 

Z=\=4- (17) 

1 + ph x ' 

The point Z=0 corresponds to ph = 1 or p'z=0, giving z = a 1 , o 2 , or o 3 , we 
may choose z = a 2 . The point Z = 1 corresponds to ph = 0, we may choose 
z = o 2 + -g- 6) 2 . The point ^"= — 1 corresponds to £> 8 s = », and we may choose 

2=0. 

The function Z = (1 — p 3 z)J(l + p s z) is finite, continuous, and one-valued 
within the right-angled triangle in the z plane whose corners have coordinates 
(0 , o 2 ) , ( o> 2 , — g w 2 J , (0,0), and the points on the real axis in the Z plane 
correspond to the points on the sides of the triangle, the points Z = , 1 , — 1 
being those which correspond to the corners. The side opposite to the angle of 
60° is of length equal to the real half-period w 2 . 



^"1 

T ° 



x 



-1 



-X 



z 



The figures are drawn as in Case I. 
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9. Suppose there is a vortex of strength mat z , and let Z be the corres- 
ponding point in the Z plane, and z' , Z' the conjugate imaginaries to z , Z , 
then the fluid motion is given by 

,,, , A wi Z—Zq m p 3 z — p\ /1oN 

* + t $ = — - log ^ » == __l g£L »L_o . (18) 

27E D Z — Zjj 27t 6 |Tz — £>% v ; 

The roots of p*z = p 3 z are z , e x z , e 3 z , say z , z lt z 2 , and those of %>h = v\ 
are the conjugate imaginaries, say z' , z{, z' % , resolving the numerator and denomi- 
nator of (18) into factors and remembering the fundamental equation (10), we 
find 

?+,*>- m l na - g ( g ~ gp)6(z— zQgfe — Z t )6(z+Z )<S(z+Z l )6{z + Z l ) / lq x 
^ w ~ 2n l0g tiz-ztt 6(z-z{)6(z- zl)6(z + z>)6(z + z[)6(z +zi)- Uyj 

Case III. — The isosceles triangle containing an angle of 120°. 

10. In this case the equation (A) becomes 

Gl z + a % =f \ z _ a)Kz d l hf{Z _ 4 ■ (20) 

Take a = 0,&=l,c = — 1, and suppose 3 = — , then 

/*' dZ' 
— (l_Z«)f 

Write , 1 £ 3 



1— ~Z n ~ Z*— 1' ( 21 ) 

1 z+G^-ZJ vw ^ Ty (22) 

Hence, adjusting constants, 

x = jpz, (23) 

where the invariants and periods of the elliptic functions are the same as in 
Case I. 

The correlation is effected by means of the equation 

Z* = - p * z 



Y*— l* 



^=v- < 24) 

The point 2T = corresponds to pz = and we may choose z= a z -\ oi . 

The points Z= ± 1 correspond to £>2 = 00. and we may choose z = for Z~= 1 
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and z = 2w 2 for Z = — 1 . Then since the function pz does not vanish anywhere 

within the triangle whose corners are z = a% + — al, 0, 2g) 2 , it follows that the 

function ph has no branch-point within the triangle, and is consequently finite, 
continuous, and one-valued within the triangle when its value has been chosen 
for one point. 

Thus, choosing one branch of the function plz, we may say that the func- 
tion Z=- 2ph/p'z is finite, continuous, and one-valued within the isosceles 
triangle in the z plane whose corners have coordinates 



(oi> 4gO|). (0. 0). ( 2 %. °)> 



■v/3 

and the points on the real axis in the Z plane correspond to the points on the 
sides of the triangle, the points Z=0, 1, — 1 being those which correspond to 
the corners. The length of the base of the triangle is the real period 2^. 

y 



i 

z 




-x 



-1 



o 
Z 



~X 



The figures are drawn as in Case I. 

11. Suppose there is a vortex of strength m at %, and let z' be the conjugate 
imaginary, Z Q , Zl the corresponding points in the Z plane. Then the fluid 
motion is given by 

y + l ®-^ [og z=z>- 2^ log lW^ Tz-J/KVz- wJ) ( } 

The function ^—. ^-2 is not a rational function of pz, p'z and does not 

p'z f'z 

appear to be expressible as a product of 6 functions. 
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Case IV. — The isosceles right-angled triangle. 

12. In this case the equation (A) becomes 

dZ 



C* + G 2 -f (z _ a)HZ -b)i(Z-c)i ■ (26) 



1 



Take a = 0, b= 1 , c = — 1 and suppose Z= -=-. , then 

Z 

dZ< 



GlZ+0 »~-f Z (i-z>*)i 



then 00+0,=- r. *> rr (28) 

J v x(ar— 1) 
Hence, adjusting constants, 

x = pz, (29) 

where $>z is defined by the equation 

p\ — 4(f% — e x ){<pz — e % )(pz — e 3 ), ) > ZQ . 

e x = 1, ^ = 0, e 3 = — 1. ) 

The correlation is effected by means of the equation 

Z* 



Z % —\ 



re, 



giving z _ $>g _ pz — e 2 ^ .^ 

Vj? 2 2— 1 V(pz — e x )(pz — e 3 ) ' 

15. In the elliptic functions given by (30) we have to notice that the 
discriminant is positive, and the invariant p s = 0. The real half-period a is 



•=JT 



\/4x(x 8 — 1) 4 p /_3 



(i) 



and the pure imaginary half-period 6>' = ia. 

One root of ^>2 is easily found since f (a + w') = e 2 = 0. 
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Again Vpz — e 1 and v^z — e 3 are singly determinate functions, viz. 



Vpz — e 1 — 



6& 



6»z 



pz 



&z 



6z 

Hence 



and \/<pz — e s = ^- , and since <s 2 = we have in like manner 



Z = 



<o 1 z<o s z 



(32) 



14. The point Z=0 corresponds to pz = and we may choose z = a + «'• 
The points Z= ± 1 correspond to pz = «, and we may choose z = for Z = 1 , 
and a = 2o> for ^= — 1 . 

The function Z= S|z/(6 x z6 3 z) is finite, continuous, and one-valued within the 
isosceles right-angled triangle in the z plane whose corners have coordinates 
(o, a), (0, 0), (2o, 0), and points on the real axis in the Z plane correspond to 
points on the sides of the triangle, the points Z — 0, 1 , — 1 being those which 
correspond to the corners. The length of the hypothenuse is the real period 26). 

y 




The figures are drawn as in Case I. 

15. Suppose there is a vortex of strength m at the point z , and let Z be 
the corresponding point in the Z plane, Z(, the image of Z in the real axis and 
zl the corresponding point to Z[, it will be necessary to prove that z , z[ are 
conjugate imaginaries. We shall show that if z' be-the conjugate imaginary 
to %, then Zq is conjugate to Z . 

Since yf and a are real, we see that the function <o x (a + i(3) , being equal to 

<o (a + iB + o) , , „ . , 

gr-iia-t-w 1S foQ conjugate imaginary to 6 X (a — iB). 

* ? is the constant -— which = \-f ( — — y>x ) dx . The constant >/ is — , and the relations 

,<y _ ,' w — _ OT , «' = ju , r/ = — m have place (Halphen, ch. V) . 
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The function 6 2 (a + t/2) = ^L + ^+ii^Ml <r ,(.+«+.,(-+* and the 

^ V ' ^' 6(&> + tO) 

exponential factor is e -i(«+«-»>(e-«). Also 

<S 2 (a - </?) = gK«+J»"'0g + ") + ^l L <-,(.-»+.,* + ., 
•* v r y 6(g) — to + 2o ) 

— 6 U<* + Q ) — ' (ft + Q )} ^V(— ^Q + ^tf— ) 

6 (a — to) 

_ 6|(a + 6)) — t(/? + o)} ^_,(. + p ) + t ,o-.) i 
6(g) — iu) 

Thus 6 2 (a — t/?) and 6 2 (a + i(3) are conjugate imaginaries. 

Again, the function 6 3 (a + t/3) = gjfL+i(ft±iii fl «(-+* and the expo- 
nential factor is <r * <" - la) . Also 



^(— w= ^i^^ M ^^ 



6 (— to + 2w') 
- 6 (— to) ^ 

_ 6|a— t(/? + o)f p -,()3 + t a) 

- 6(— to) 

Thus 6 3 (a + t/2) and 6 3 (a — t/3) are conjugate imaginaries, and hence it 
follows that z , 2o being conjugate imaginaries, Z 0) Z' are so, and inversely. 
16. The fluid motion is given by 

Kemarking that pz — 1 = <o\zj<oh , 

and fpz + 1 = 6|s/(3 a s, 

so that 6?z + <o\z — 2&zpz = 2fz , 

equation (33) becomes 

m. rh_ m 1 (glZglZp — <0 Z Z6 3 Z )(6 1 Z<0 Z Z (S — ggggigg) (o 4 \ 

* + l ^ - 2ti 10g (6 lS (3^ - cveWXMvtf - <5*<^) ' 
The factor ^z^gZo — 6 1 2 6 3 z of the numerator may be put into factors. Using 
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Schwarz's formula [D][9], p- 81, and taking first [i=l, /L=3, and then 
(i =. 3 , h = 1 , we have 



r- /- ^2 + 20^2+20^2 2 ,- 2 2n 

S1Z0S32 = <*i -^-5 6 3 — L_S 6 X — __» 63 _^ 



2 +_2o g 2+2 g 2 — 2 g 2 — 2 

2 2 2 2 a 2 



whence 

6,26320 - 6 lZo 6 3 z = 2 ( ei - e 3 ) <3 2 -±^ 6 z -=^ 6, £±* 6 2 ^p> , 

So that the fluid motion is given by 



TO 



(S^CVo— 6 3 z6 8 zo) 6 — (2+20) 6 — (2— a ) 6 2 — (2+2 ) 6 S — (2— 2 ) 



<P+ t * = ^ log 1 I 1 1 , (35) 

(6^6^— 6 8 26 3 2') 6 _ (2+2') 6 — (2— 2 ') 6 2 ~2 ( s +«o) 6 3 — (z—z'o) 

17. It may be noticed that in the notation of Jacobi's elliptic functions 

-^= ,«=-—. (35a) 

cnz \/2 v y 



So that 



¥ + t <l> = -|Llog 



/ 1 



-7r(l+cn a 2) 



\ cnz 

2 ^\ //— ri4-<m*^ 



1 (l+cn*2 )\//4- (l+cn 2 2) i- (l + cn»«0 \ 



en z // \ en 3 en 2« / 

rp + ^ = -^log { Cng ~ CnZ i ( ?~ CngCn 4, (36) 

2n (en 2 — en 2«)(1 — en 2 en 20) v y 

is the solution for a vortex at z . 

As an example of the application of the method when the motion is due to 
sources or sinks, the solution for a source and sink each of strength m at the base 
angles marked 1 and — 1 is 

* + *=«.! gi+^_ 3 . (37) 

n & 1 — en 2 
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